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ABSTRACT 

The effect of gravitational lensing on cosmic microwave background (CMB) 
anisotropics is investigated using the power spectrum approach. The lensing 
effect can be calculated in any cosmological model by specifying the evolution of 
gravitational potential. Previous work on this subject is generalized to a non-flat 
universe and to a nonlinear evolution regime. Gravitational lensing cannot 
change the gross distribution of CMB anisotropies, but it may redistribute 
the power and smooth the sharp features in the CMB power spectrum. The 
magnitude of this effect is estimated using observational constraints on the 
power spectrum of gravitational potential from galaxy and cluster surveys and 
also using the limits on correlated ellipticities in distant galaxies. For realistic 
CMB power spectra the effect on CMB multipole moments is less then a few 
percent on degree angular scales, but gradually increases towards smaller scales. 
On arcminute angular scales the acoustic oscillation peaks may be partially or 
completely smoothed out because of the gravitational lensing. 
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1. Introduction 



The effect of gravitational lensing on the cosmic microwave background anisotropies 
has been studied in the past by several groups ([Kashlinsky 1988| ; | Blanchard fc Schneider 



[1987| ; Cole fc Efstathiou 1989j [Sasaki 1989| ; [lbmita fc Watanabe 1989t |Linder 1990a| ,b 



Cayon, Martinez-Gonzalez & Sanz 1993a,b; [Bassett et al. 19U4]; Fukushige et al. 1994). 



Using different approaches these authors came to very different conclusions about the 
importance of the effect. Cole & Efstathiou (1989) used a nonlinear CDM model and 
found a small effect on CMB. Cayon et al. (1993a) used a linear model and found an 
appreciable effect on arcminute angular scales for some models. On degree angular scales 
they also found a negligible effect. On the other hand, using different approaches such as 
the Dyer-Roeder distance-redshift relation or simplified N-body simulations Bassett et al. 
(1994) and Fukushige et al. (1994) found a significant effect even on degree angular scales. 

There are several shortcomings of these studies that do not allow one to draw a firm 
conclusion on the importance of the lensing effect on CMB. First, the studies are based 
on a particular cosmological model and the results could change significantly if the model 
is changed. While some groups (e.g. |Cole fc Efstathiou 1989| ; Cayon et al. 1993a,b|) 



attempted to asses this uncertainty by presenting results for different viable cosmological 



models, others (e.g. Passett et al. 1994] , Fukushige et al. 1994) used models that do not 
allow a direct comparison with existing observational constraints and thus may not even 
be realistic models. A second shortcoming of previous studies is that they do not fully 
include the evolution of large-scale structure in their models. While Cole & Efstathiou 
(1989) calculated the effect only at late epochs when the matter is in the nonlinear regime, 
Cayon et al. (1993a,b) only included the linear evolution, whereas Bassett et al. (1994) and 
[Fukushige et al. (1994)| neglected any evolution at all and assumed that the universe did 
not change from a certain redshift until today. 

The purpose of this paper is to provide a more realistic answer on the importance of the 
effect by using observational constraints on large-scale structure distribution and properly 
including its evolution. The method used is based on the power spectrum approach in 
linearized gravity and is similar to the one used by Linder (1990a,b), Kaiser (1992) and 
Cayon et al. (1993a,b). An equivalent method based on optical scalars has been developed 
by Gunn (1967) and extended by Blandford & Jaroszynski (1981). Present work differs 
from previous studies in that I also include the nonlinear effects by modeling the power 
spectrum evolution in the nonlinear regime. By comparing the nonlinear calculation to the 
linear approximation one can identify the angular scale where the nonlinear effects become 
important. I also extend the calculation to the case of an open (or closed) universe and 
correct some erroneous expressions in the literature, all of which allows to calculate the 
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lensing effect in any standard cosmological model (i.e. in any model based on a weakly 
perturbed metric in a universe that is homogeneus and isotropic on large scales). The 
estimate of the lensing effect on the CMB is based on the observational constraints on the 
power spectrum and on the ellipticity correlations of distant galaxies, which enables to asses 
its magnitude in our universe. The results are presented in terms of the CMB anisotropy 
power spectrum, which allows one to discuss the effect independent of the observational 
strategy. In the conclusion section I discuss the possible sources of discrepancy between 
present results and some of the previous work on this subject. 



2. Formalism 



In this section I review the formalism to compute the gravitational lens effect on 
a pair of propagating photons separated by an angle 6 at the observer's position. The 
starting point is a perturbed Robertson- Walker model with small-amplitude scalar metric 
fluctuations. In the longitudinal gauge ( [Bardeen 1980| ) one can write the line element as 
(adopting c=l) 



ds z = a 2 (r 



-(1 + 2<p)dr 2 + (1 - 20) [dx 2 + sm 2 K X (d6 2 + sin 2 



where the metric is expressed with comoving spherical coordinates and conformal time r 
and a(r) denotes the expansion factor. I defined 

K-^sinK^x, K > 
sm /v x = { X , K = (2) 
[-K)-V 2 ±}h(-K)V 2 x, K<0 

Curvature K can be expressed using the present density parameter Q and the present 
Hubble parameter Hq as K = (Qq — 1)H 2 . The total density O(r) is in general time 
dependent and can have contributions from mass density fi m (r) or vacuum energy density 
Q v (t), fi(r) = f) TO (r) + Q v (t). The metric perturbation <f> can be interpreted as the 
Newtonian potential since, neglecting the contributions from wavelengths larger than the 
Hubble distance, it obeys the cosmological Poisson equation 

V 2 = ^H 2 n m (r)a 2 6, (3) 

where S is the mass density fluctuation. Statistical properties of the potential on scales 
small compared to the curvature scale can be described with its Fourier transform 0(fc,r), 
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where <p(x, r) = / d 3 k<f)(k,T)e tkx . Its ensemble mean and variance are (4>(k,r)) = and 

r)(p*(k f , r)) = P<^(A;, r)<5 3 (A; — A;'), where P^(A;, r) is the power spectrum of the potential 
at time r. 

A photon propagating through the universe will be deflected by the mass concentrations 
along its path. The rate of change in the photon direction n is given by the photon geodesic 
equation, which applied to the metric in equation [l] gives 

dfl -* -* 

— = 2n x (n x V0) = -2V ± <j>, (4) 
al 

where the symbol Vj_0 denotes the transverse derivative of the potential and I is the 
comoving path length along the photon geodesic. Gravitational potential <fi can be viewed 
as providing a force deflecting the photons while they propagate through the unperturbed 
space-time, described by a 3-sphere (closed universe), 3-hyperboloid (open universe) or 
Euclidean space (flat universe). Because the only observable photon direction is that at the 
observer's position it is convenient to propagate photons relative to their final direction (i.e. 
backwards in time). Gravitational lensing is not expected to lead to large deflection angles 
(e.g. Linder 1990; Seljak 1994) and one can replace the transverse derivatives in equation 
(|j) with the transverse derivatives with respect to the observed direction of the photon or 
with respect to any other direction that has a small angular separation with the photon. In 
this spherical plane approximation the observed direction of the photon can be described 
with a two-dimensional angle 9 with respect to the origin. Moreover, the null geodesic 
condition for photons gives dr rs — d\ neglecting corrections of the order 0(0). Even when 
metric perturbations are present, one can continue to parametrize the geodesic with the 
unperturbed comoving radial distance \ or the conformal time r, which are related through 
X = t — r, where r is the conformal time today. 

The total deflection angle between the photon source at the last-scattering surface 1 
and the observer is given by 

5a = -2 [ XlCC V ± <pdx- (5) 
J o 

Similarly, the photon angular excursion on the last-scattering surface relative to its observed 
value is given by 

-. /"Xrec _, 

59 = -2 W(x,Xr CC )V±Hx, 
Jo 



1 I assume throughout the paper that CMB anisotropics are generated at a recombination time r rec ~ 
and therefore x roc ~ tq. This is a good approximation for all the models where the CMB fluctuations are 
generated at a high redshift. 
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W{X,Xrec) = : • (6) 

Sm K Xrec 

In a flat universe the latter simplifies to W(x, Xrec) = (1 — x/Xrec)- Note that it is 88 that 
is relevant for the discussion of lensing effects on CMB, because one is interested in the 
angular excursion of a photon on the CMB last-scattering surface and not in the change in 
its direction. Some of the previous work on this subject used 8a instead of 88 (e.g. |Cay6n 



st al. 1993a, b|). As shown by Seljak (1994), in a flat Q m o = 1 linear theory this leads to a 



factor of 10 1 / 2 overestimate of the relative dispersion between two photons. In the following 
I will restrict the discussion to 86. 

Two photons A and B observed with an angular separation 6 have a different angular 
separation when emitted from the source position. Its mean is equal to the unperturbed 
value, while the dispersion is given by (Seljak 1994), 

a{6) = 2- 1 / 2 (^86 A -86 B y /2 = C gl (0) - C gl {6) 

POO /"Xrec 

C gl (6) = 167T 2 / k 3 dk P^k,T = T - X )W 2 (XiXrec)Jo(kdsm KX )dXi (7) 
Jo Jo 

where ()# denotes the ensemble average performed over all pairs of photons with a fixed 
observed angular separation 6 and Jq(x) is the Bessel function of order 0. The derivation 
of equation [1] is based on Limber's equation in Fourier space (e.g. [Kaiser 1992j ), which 



assumes that the dominant scales contributing to the dispersion are much smaller than the 
photon travel distance. This condition is satisfied for sources at cosmological distances. No 
assumption on the power spectrum has been made and equation ^ can be used both in the 
linear and in the non-linear regime, both of which can be described by a time evolution 
of the power spectrum P^(fe, r). I neglected the Poisson contribution arising from the 
discrete nature of galaxies, which is only important if Q m o < 0.1 ( piandford fc Jaroszynsk]| 
1981 ; Cole fc Efstathiou 1989|) . An important approximation needed to derive equation [?| 



is that the potential sampled by the perturbed photon geodesic can be replaced with the 
potential along the unperturbed geodesic. This limits its applicability to the regime where 
&(&)/& <C 1. When u(6)/6 becomes of order unity the above assumption is no longer valid 
and the separation between the photons may grow significantly faster than what the model 
would predict. This point will be discussed again in the last section. 

It is useful to give a physical understanding of the lensing effect on two nearby photons. 
For simplicity I restrict the discussion to scattering on a single scale k~ x and to a flat 
space. In real space k~ x is a correlation length and determines the scale on which regions 
become uncorrelated. For sufficiently small angles the photons separated by distance x@ 
propagating through a region of size /c _1 are coherently deflected (assuming x6 < k~ l ). 
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The change in the relative angle between the two after crossing this region is given by 
59% ~ X^V^(2A; _1 Vj_0) ~ 2kx<p9. The photons pass through iV m kxrec uncorrelated 
regions and the total rms deflection angle between the photons grows in a random walk 
fashion, 59 = N l l 2 59i or a{9)/9 oc (/cx rec ) 3//2 0. Adding the contributions from different 
modes one reproduces, numerical factors aside, the small angle limit of equation 0. For 
large separation angles the scattering is incoherent (xO > k~ x ) and each photon is deflected 
by (^Xrcc) 1 / 2 ^ independently of 9, implying a{6)/6 — > 0. This asymptotic behavior is 
confirmed by the numerical results presented in the next section. 

Once a(9) is known as a function of 9 it is straightforward to calculate the lensing effect 
on the CMB fluctuations. The effect is most easily expressed in terms of the temperature 
anisotropy correlation function C{9) — (^sr 4?- )$■ Using the two-dimensional formalism 
and isotropic approximation presented in the Appendix one obtains the modified correlation 
function C{9), 

9(3 



a 2 (9) Jo 



a 2 (9) 



where Jo is the modified Bessel function of order 0. This equation is strictly valid only for 
gaussian fluctuations, but should give a reasonable estimate of the effect even when this 
condition is not satisfied. Note that the effect of lensing is to integrate the correlation 
function with approximately a gaussian centered at 9 with dispersion <j(9), as can be seen 
using the asymptotic expansion of Iq combined with the exponential in equation §. Thus, 
lensing acts as a filter smoothing out the sharp features in the correlation function. For 
lensing to be important the correlation function at 9 must be changing rapidly on a scale 
a{9). 

It is customary to present different models of CMB anisotropies in terms of the 
power spectrum, given by the multipole moments C\. These are obtained by the Legendre 
transform of correlation function, C\ = 2n Jq sin 9C(9)Pi (cos 9)d9, where Pi(x) is the 
Legendre polynomial of order I. The lensing effect on the G\ multipoles can be efficiently 
calculated using the Gauss-Legendre integration of C{9) in equation [X6[ One can estimate 
the effect on C\ by assuming e = cr{9)/9 is a constant (Bond 1995). While this is not true 
in general (figure [ft), one may hope to use this approximation if a{9)/9 is sufficiently slowly 
changing with 9 and is small, so that only correlations over a narrow range of 9 are mixed 
by lensing. In this case e should be determined by its value at 9 oc / _1 ( [Bond 1995[ ). From 
equation [AO follows in the e<l limit 

q= fc ( ,-l e -H')V^f. (9) 



Lensing thus smoothes the spectrum of Ci with a gaussian of relative width e, similar to the 
effect on the correlation function. 
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One way to calculate a(9) is to use the observational constraints on the power spectrum 
from the large-scale structure observations, carefully including the effects of the evolution 
in a given cosmological model. This approach will be explored in the next section. A 
somewhat less model-dependent estimate can be obtained from the observational constraints 
on correlated distortions of distant galaxy images. This can be described by p(9), the 
average polarization within a circular aperture of radius 9, which describes the correlations 
in the ellipticities of galaxy images as a function of angle 9. It is related to the power 
spectrum using a Q ^ 1 generalization of the expression given by Blandford et al. (1991) 
and Kaiser (1992), 

p 2 (9) = lQ7r 2 k 5 dk P4k,r = r - X ]W 2 (x,X g )sm 2 K x 
Jo Jo 

where I assumed for simplicity that all the galaxies lie at the same source position Xg- Using 
a small argument Taylor expansion of Bessel functions in equations |7| and [H] one obtains 
a simple scaling between a (9) and p(9) for a flat Q m o = 1 universe in the linear regime, 
independent of the power spectrum on small angular scales, cr{9)/9 = 2~ l p{9) (x rec / 'x g ) 3 ^ 2 ■ 
Nonlinear effects and Q m o < 1 make the low redshift contributions more important relative 
to the case above, which decreases a(9)/9 derived from p(9). Numerical evaluation confirms 
this prediction and so the scaling above can be used to give an upper limit on a(9)/9 from 
the observational limits on p(9) on arcminute scales. 



2Ji{k9 sin^ x) 
k9 sin^ x 



d X , (10) 



3. Estimate of the Lensing Effect in Our Universe 



To compute the lensing effect one needs to specify the power spectrum of potential as 
a function of scale and time. In linear regime the time dependence of density perturbation 
in a CDM dominated universe obeys the well known growing mode solution. For the 
particular case of Q m o = 1 universe the potential does not change in time and lensing 
contributions at early times are as important as those at late times. For the nonlinear 
evolution of the power spectrum one can either use N-body simulation results or adopt a 
semianalytic approximation for it. Here I adopted the approximation given by [Hamilton] 
|et al. (1991)| , generalized to Q ^ 1 by Peacock fc Dodds (1994)1 and to the density power 
spectra with slopes n < — 1 by [Mo, Jain fc White ( 1995 )| . This prescription is based on an 
educated guess of what the evolution of the density correlation function should be in the 
nonlinear regime. Although not exact, it agrees well with the results of N-body simulations 
(see Peacock & Dodds 1994 and Mo et al. 1995 for a detailed discussion of its applicability) 
and should give a good estimate of the nonlinear power spectrum in the regime where 
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dissipative baryonic processes can be neglected. The linear to nonlinear mapping is most 
easily expressed using the mass density variance A 2 (A;), which is related to the potential 
power spectrum via A 2 (k) = 167rfc 7 P^(A;)/9fi m (r) 2 iJ 
and nonlinear power spectrum is given by 



4 a 4 . The relation between the linear 



A 2 (k nl ) 

fnl(x) 



f nl [A 2 (kt)]; k l = [l + A 2 (k nl )}- 1 / 3 k nl 
1 + 0.2px + (Ax) a/3 



x 



[1 + ([Ax] a g 3 (n m , Q v , a)/[11.68xV2])/3 



(11) 



where A = 0.84[g(fi m , Q v , a)] 0,2 , a = 2/g(fl m ,Q v ,a) and (3 = 2g(Q m ,Q v ,a). The linear 



growth factor g(Q m ,Q v ,a) can be approximated with a few percent accuracy as ( [Lahav eT 
al. 1991| ; |Carroll, Press fc 'l'urner T99p 



g(Q 



X 



5Q 



rnO 



2[X(1 + [{l m0 /aX]™) - a 2 n v0 + n m0 /2a] 
1 + n rn0 (a- 1 - l) + Q v0 ( 



(a 2 - 1) 



(12) 



Mapping in equation |ll] can be improved by allowing for the variation in the shape of 
the power spectrum. One can introduce an effective index of the density power spectrum 
n cS = c?ln[A; 4 P^(/c)]/(iln/c at a wavevector k defined such that the rms mass fluctuations 
averaged within a sphere of size k" 1 is unity (cr^-i = 1). For n c g = the mapping above 
gives reliable results, while for n e R < — 1 there are substantial deviations from the N-body 
simulations. As shown by Mo et al. (1995), one can improve this by replacing A 2 by 
A 2 / B(n ef [), where 



B(n) = 0.795 



r 


' 17+n \ 


K 10+2n J 


r 


( 11+n \ 


Vl0+2n J 



-(5+n) 



(13) 



In this case a better fitting formula for the Q m o = 1 case is given by (Mo et al. 1995) 

/l + 2x 2 -0.6x 3 - 1.5x 7 / 2 + lx 4 \ 1/2 



fm(x) 



X 



1 + 0.0037x 3 



(14) 



The system of equations presented above can be used to calculate u(9)/9 for most 
cosmological models of current interest (one exception being the models with massive 
neutrinos on small scales where neutrino free streaming is important). To obtain an 
estimate of the lensing effect in our universe I will use observational constraints on the 
power spectrum, as compiled by |Peacock fc Dodds (1994)| . The power spectrum can be 
parametrized with a CDM type linear transfer function ( Pardeen et al. 1986Q with two free 
parameters, the amplitude a 8 , determined by the mass fluctuation averaged within a sphere 
of radius 8/t _1 Mpc and the shape parameter Q m oh, determined by the turnover position 
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in the power spectrum (h is the present day Hubble parameter in units of lOOkm/s/Mpc). 
For wavevectors between 10~ 2 and 1 /iMpc -1 all the galaxy and cluster surveys are in a 
reasonable agreement with a CDM type linear power spectrum with Q m oh 0.25 (Peacock 
& Dodds 1994; da Costa et al. 1994 ). For normalization I will adopt <jg = 0.8, which is close 



to the normalization obtained by [Peacock fc Dodds (1994)] and by [White, Efstathiou fc 



Frenk (1993)| using the cluster abundances normalization over most of the interesting range 
of flmo- This normalization also agrees with the COBE normalization for the currentlty 
favored Q m o = 0.4, h = 0.65 case, both in the open universe model ( porski et al. 1995| ) 
and in the cosmological constant dominated model with a modest tilt (Stompor, Gorski & 
Banday 1995). The adopted power spectrum is likely to be within a factor of two of the 
real power spectrum on the arcminute scales and larger. 

Given the linear power spectrum and its nonlinear evolution one can compute 
e = a(8)/8 as a function of 9. In figure [l] the results are presented for the power spectrum 
discussed above in flat Q m o = 1, flat low fl m o model and open low Q m o model, all of which 
are known to phenomenologically agree with most of the large-scale structure observations. 
The thick curves give the result of a full nonlinear calculation, while the thin curves show 
the corresponding linear case. One can see that while in the linear case a{9)/9 approaches 
to a constant for small 9, it continues to increase in the nonlinear case. Therefore in the real 
universe one cannot define a typical coherence angle, which was used by previous studies 
(e.g. |Sasaki 1989| ; [Linder 1990a| ; |Uay6n et al. 1993a, b| ) and the approximation e = const 



is not valid on any scale. The results can only be used on angular scales above a few 
arcseconds, where a(9)/9 <C 1 and where the nonlinear mapping (based on the evolution of 
collisionless matter) gives reliable estimates. 

As seen in figure [L] in the linear regime the lensing effect on CMB decreases with fl m Q. 
This is mainly due to the linear decrease of potential with Vt m in the Poisson equation, partly 
offset by the longer travel distance, larger growth factor ratio g(fl m , &v, a )/g(^mo, &vo, a o) 
and in the nonlinear regime by larger nonlinear effects in the low Q m o models. The latter 
is more important because the scales that are nonlinear today became nonlinear earlier 
than the corresponding scales in an Q m0 = 1 universe. In a low Q m Q model the universe 
changes from Q m Q « 1 to Q m o <C 1 earlier than in a cosmological constant model with the 
same matter density and in addition the relation between the conformal time and angular 
distance changes, all of which leads to a larger lensing effect on very small scales. The value 
of a (9)/ 9 linearly increases with <jg in the linear regime, but grows faster than that in the 
nonlinear regime. While the value of <7g is still somewhat uncertain, it is unlikely that erg is 
much bigger than 1 even in an open model and the curves on figure [I] should indicate the 
range of the lensing effect in our universe. To investigate the sensitivity of the effect to the 
shape of the power spectrum I compared the flat model above to the standard CDM model 
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Fig. 1. — cr(9)/9 versus 9 for 3 different values of f2 m0 and fl v0 using the power spectrum 
with Q m o = 0.25 and ag = 0.8. Thick lines are the result of a full nonlinear calculation, while 
the thin lines give the corresponding linear case. Also indicated are the 90% c.l. upper limits 
from ellipticity correlations of distant galaxies, as derived from observations by Fahlman et 
al. (1994) (A) and Mould et al. (1994) (B). 



with Q m0 = 0.5. The relative difference between the two models at 9 ~ V only depends on 
the power spectrum amplitude, because the dominant scales there are similar to the scales 
that contribute to the <jg normalization. The inverse wavenumber that makes a dominant 
contribution to a(9) is approximately 0.5 h Mpc for 9=1' and 0.05 h Mpc at 9 = 1° and 
above. On larger angular scales the differences in a(9)/9 between the different spectral 
shapes increase (with standard CDM model having less power and thus smaller a(9)/9 for 
a given ag normalization), but the overall effect is decreasing and becomes rather small on 
degree angular scales independent of the model. 

In figure p] the 90 % c.l. upper limits on u(9)/9 are also indicated, as derived from 
Mould et al. (1994)| and [Fahlman et al. (1994)| limits on the correlated ellipticities. Both 
groups report a null detection of average ellipticity within a 4.8' and 2.76' radius aperture, 
respectively, with a sensitivity of about 1%. Adopting median redshifts of z = 0.9 and 
z = 0.7 gives radial distances 0.27 and 0.23 times the comoving distance to the horizon, 
respectively. For the two surveys one obtains upper limits that are comparable to the power 
spectrum estimates, which gives additional confidence that the effect was not severely 
underestimated. A general conclusion that can be derived from these results is that cr(9)/9 
is less than 20% on scales above 1' and less than 5% on scales larger than 1°. 
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Fig. 2. — CMB anisotropy power spectrum 1(1 + l)Ci versus I with lensing (dashed lines) 
and without lensing (solid lines). Upper curves are for adiabatic CDM model with h = 0.5, 
O m0 = 0.4 and Q v q = 0.6, lower curves are for adiabatic CDM model with h = 0.5, fl m o = 1 
and Q v o = 0. Both models are normalized to COBE. Lensing smoothes the sharp features in 
the power spectrum, but leaves the overall shape unchanged. The two models show a typical 
range of the lensing effect on CMB. 

Figure |2] shows the lensing effect on the CMB fluctuation power spectrum for the 
models discussed in the previous paragraph. The CMB multipole moments were obtained 
from a numerical integration of perturbed Einstein, Boltzmann and fluid equations ( [Bode 
|fc Bertschinger 1995| ). These models exhibit characteristic acoustic oscillations (Doppler 
peaks) and suppression on small scale due to the diffusion (Silk) damping. Lensing 
induces very little gross change in the power spectrum of CMB. However, the peaks of 
acoustic oscillations are smoothed because of lensing and on smaller angular scales they 
can be completely erased. This occurs both because cr(0)/9 increases and because the 
relative width of the oscillations becomes narrower towards the smaller angular scales. 
Observational sensitivity to this effect depends on the particular experimental setup, but 
for most experiments the window functions are relatively broad in /-space and consequently 
the effect is diluted. 

Currently popular models predict that primary CMB anisotropics are only important 
above the Silk damping scale of order of a few arcminutes. As shown in figure [1], on these 
scales the linear contributions to the dispersion <j(9)/9 are still dominant and nonlinear 
effects are negligible. In fact, using only the linear theory evolution of power spectrum 
gives results indistinguishable from the fully nonlinear calculation in the regime of interest 
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Fig. 3. — Comparison between several approximations for calculating the lensing effect on the 
CMB anisotropies in the COBE normalized CDM model with h = 0.5 and fl m oh — 0-5- Both 
nonlinear and linear isotropic approximations gives results that are almost indistinguishable 
from the fully nonlinear and nonisotropic calculation over this angular range, while e = const 
approximation gives reliable results only over a limited range of I and cannot be used for an 
accurate calculation of the lensing effect. 

(I < 2000, see figure |[). These angular scales are thus unaffected by the uncertainties of 
the nonlinear evolution and are also the scales where the large-scale observations place 
the best constraints on the power spectrum. The approximation given in equation |9| gives 
reasonable results only if one adopts e at 9 « AttI^ 1 , which is significantly larger angle than 
expected (Bond 1995) and again implies that nonlinear effects are not important until very 
large I. Even then the agreement is only approximate (figure |3]) and limited to / < 1000. 
In general it is better to use the isotropic approximation in equation together with the 
Gauss- Legendre integration to calculate the lensing effect on the multipole moments, as it 
is not significantly harder to compute than the approximation given in equation |9|. 



4. Discussion 



The lensing effect on the primary CMB anisotropies can be calculated for any 
cosmological model with a specified evolution of gravitational potential power spectrum. 
This formalism was applied to several currently popular models, which best fits the 
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observational data on large scales. The results indicate that the gravitational lensing does 
not significantly affect the CMB power spectrum on degree scales and larger, but becomes 
gradually more important towards the smaller scales. Lensing redistributes the power in 
the angular correlation function and the amplitude of the effect depends on the smoothness 
of the underlying CMB spectrum. For standard adiabatic models the acoustic oscillation 
peaks are rather prominent even at small angular scales (beyond / ~ 1000) and lensing may 
completely erase this structure. 

Recently, two groups claimed that the gravitational lensing effect on CMB has been 
severly underestimated in previous calculations and that it importantly changes the CMB 
pattern even on degree angular scales. Bassett et al. (1994) assume a model in which 
photons propagate through a homogeneous universe with a density smaller than its mean 
density to account for the fact that some of the mass resides in dense clumps. Using 
the Dyer-Roeder distance-redshift relation in such a universe they obtain an increase in 
angular separation between the two photons relative to its unperturbed value. Similarly, 
Fukushige et al. (1994)| assume a model in which the universe is populated by a number of 



massive clumps embedded in a large empty void. Here the angular separation between two 
photons is additionally increased with every passage of the photons by a clump, because 
the closer photon is always deflected more than the more distant one. This leads to an 
exponential growth of the angular separation until it reaches the mean projected separation 
between the clumps. One problem with these models is that they cannot be applied to the 
large (supercluster) scales, where the density fluctuations are small, given that in these 
models the mass density in a box on scales smaller than the mean distance between the 
clumps is either zero or very large. Observational data on large scales suggest that density 
fluctuations are close to gaussian and both underdensities and overdensities have to be 
included for a proper description of light propagation. The effect of underdensities is to 
decrease the angular separation between the two photons and this leads to a random walk 
growth of rms deviation between them. As shown in this paper the lensing effect on CMB 
on arcminute scales and larger is dominated by the linear regime, where underdensities 
and overdensities play equivalent roles. Numerical studies of light propagation in realistic 



models ( Jaroszynski et al. 1991 ) confirm that there are no large distortions in the relative 
photon trajectories present for most lines of sight, at least on scales above their resolution 
scale of a few arcseconds. Another problem with the above models is that flux conservation 
requires that the angular distance between the photons (as defined in a homogeneus 
universe with the same density parameter) remains on average unchanged, implying that 
exponential growth in separation between photons passing on the same side of a clump is 
balanced by a strong focusing of photons that pass on the opposite sides of the clump. This 
would lead to multiple images (strong lensing), but it is known observationally that such 
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situations are rare in our universe, especially on very large scales, such as superclusters. It 
is nevertheless possible that in a highly nonlinear regime (i.e. on very small scales) our 
universe could be well approximated by the models discussed by |Fukushige et al. (1994)| and 
Bassett et al. (1994). In such a regime the models presented in this paper should predict 
relative fluctuations larger than unity and their predictions would become unreliable, 
because the assumption that the potential deflecting the photons can be calculated along 
the unperturbed paths would not be satisfied. As long as u(9)/9 remains small this is 
not the case and for the density fluctuations as measured in our universe this condition is 
satisfied at least on angular scales above a few arcseconds. 

Although gravitational lensing is of small significance for the present day experiments, 
mostly sensitive to the degree angular scales, it may become relevant for the future 
experiments that will probe smaller angular scales with a much higher sensitivity and sky 
coverage. Gravitational lensing effect will be especially important for the high precision 
determination of cosmological parameters planned for the next generation of experiments. 
The uncertainties caused by the gravitational lensing should be included in the modelling 
of extraction of cosmological parameters from the CMB measurements. The formalism 
developed in this paper allows to calculate the lensing effect on the CMB for any specified 
cosmological model and can be included as a postprocessor to the standard calculations of 
the CMB multipole moments. 

I would like to thank Ed Bertschinger, Dick Bond and Enrique Martinez-Gonzalez for 
useful discussions. This work was supported by grant NASA NAG5-2816. 



A. Appendix 

Calculating the gravitational lensing effect on the CMB can be cumbersome in 
general, but it simplifies considerably if only small angular scales are considered and if the 
fluctuations in relative separation between the two photons can be considered gaussian. The 
first assumption is not very restrictive, since one does not expect the lensing effect to be 
important on large angular scales. The second assumption should really limit the validity of 
the calculation to the linear scales only, where the prediction of most models that the initial 
fluctuations are gaussian guarantees its validity. In reality its validity extends beyond that 
to the quasi-linear scales, because the relative fluctuations are obtained by a projection of 
a 3-dimensional distribution over a broad radial window function and are in general much 
more gaussian than the 3-d distribution of the gravitational potential derivative itself. 
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In the spherical plane approximation one can write the temperature anistropies ^-(9) 
in terms of its Fourier transform, 
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The correlation function including lensing is given by 
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where 59 and 59 are the angular excursions of the two photons that are separated at the 
observer's position by an angle 9. In equation |A2| one has to average both over the intrinsic 
temperature anisotropies ^r(i) and over the lensing fluctuations 59. The first averaging 
gives the angular power spectrum of CMB, 
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while the second gives the characteristic function of a gaussian field 59 A — 59 B , 
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Performing the ensemble averaging over the lensing fluctuations 59 in the equation above 
with the help of Limber's equation (Kaiser 1992) leads to the lensed correlation function 
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where a (9) is the rms dispersion in the angular positions of the two photons defined in 
equation |7| and C g \^{9) can be obtained from C g \{9) in equation ^ by replacing J with J 2 . 
Assuming l 2 C g \^{9) 1 (which follows from assuming a(9)/9 1, together with I ~ 9~ l 
and C g \^{9) < cr 2 (9)) one may Taylor expand the exponential in equation [A5] and integrate 
over ipi. This leads to 
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One approximation often used in the literature is to keep only the isotropic term 
Jq(W) in equation [A6| , which gives the dominant contribution to the lensing effect. In the 
following I will assume this approximation, because it gives a good agreement with the 
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exact calculation, as shown in figure |3|. With this approximation the lensed correlation 
function becomes 

poo 

C(e) = (2n)- 1 Idle-^^^QJoiW). (A7) 



o 



This equation is essentially the same as the Wilson & Silk (1981) expression for the 
correlation function observed with an instrument that has a gaussian beam profile, the 
only difference being that in the present case the dispersion a{9) depends on the angular 
separation 9. After another Fourier transform and using equation 6.615 from Gradshteyn & 
Ryzhik (1965) one obtains equation Alternatively, one can also express the lensing effect 
directly in terms of the CMB power spectrum, 

Ci= 9d9 l'dl'e- a ^ l '*/ 2 C v JM9)JM'9). (A8) 

JO JO 

The above expression can be further simplified if one assumes that e = a (9)/ 9 is a 
constant. Equation |A8| then reduces to 
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where Iq{x) is the modified Bessel function of order 0. Further assuming e < 1 one can 
asymptotically expand the modified Bessel function Iq{x) and use V ~ / everywhere except 
in the exponential. This finally leads to equation |9|, which is similar to the expression given 
by Bond (1995). 
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